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ABSTRACT 

In 1988 Bardeen has suggested a pragmatic formulation of cosmological perturbation theory which is 
powerful in practice to employ various fundamental gauge conditions easily depending on the character 
of the problem. The perturbation equations are presented without fixing the temporal gauge condition 
and are arranged so that one can easily impose fundamental gauge conditions by simply setting one of 
the perturbation variables in the equations equal to zero. In this way one can use the gauge degrees of 
freedom as an advantage in handling problems. Except for the synchronous gauge condition, all the 
other fundamental gauge conditions completely fix the gauge mode, and consequently, each variable 
in such a gauge has a unique gauge invariant counterpart, so that we can identify the variable as 
the gauge-invariant one. Here, we extend Bardeen's linear formulation to fully nonlinear order in 
perturbations, with the gauge advantage kept intact. Derived equations are exact, and from these we 
can easily expand to higher order perturbations in a gauge-ready form. We consider scalar- and vector- 
type perturbations of an ideal fluid in a flat background; we also present the multiple components 
of ideal fluid case. As applications we present fully nonlinear density and velocity perturbation 
equations in Einstein's gravity in the zero-pressure medium, vorticity generation from pure scalar- 
type perturbation, and fluid formulation of a minimally coupled scalar field, all in the comoving 
gauge. We also present the equation of gravitational waves generated from pure scalar- and vector- 
type perturbations. 

Subject headings: cosmology: theory — large scale structure of universe 



1. INTRODUCTION 

The Friedmann world model, based on assuming spa- 
tial homogeneity and isotropy in Einstein's gravity, is 
widely accepted as a successful cosmological model, en- 
during 90 years of observational and theoretical advances 
after Friedmann's original proposition in 1922. Main 
observational and theoretical advances have been made 
in the angular anisotropies of cosmic microwave back- 
ground radiation and in the position and motion of large 
scale galaxy distribution. Relativistic perturbation the- 
ory is important in providing crucial testing grounds for 
matching the theoretical predictions with the observa- 
tion, and for the theoretical explanation of the observed 
phenomena. The relativistic linear perturbation theory 
and the Newtonian exact and nonlinear perturbation the- 
ory are generally accepted to be successful in the current 
paradigm of modern physical cosmology. This work is 
concerned with relativistic fully nonlinear and exact per- 
turbation formulation in the Friedmann world model. 

The cosmological linear perturbation theory in the 
Friedmann's world model was pioneered by Lifshitz in 
1946. Lifshitz's analysis was made in a certain gauge 
(coordinate) condition known as the synchronous gauge 
(Landau & Lifshitz 1975). A disadvantage of the syn- 
chronous gauge condition was that even after imposing 
the gauge condition there remains remnant gauge (co- 
ordinate) mode which needs to be traced carefully. In 
this way the algebras become unnecessarily complicated. 
Often, removing the remnant gauge mode in the syn- 
chronous gauge merely corresponds to taking just an- 
other gauge condition; for example, setting the velocity 
component of the pressureless matter (being a gauge- 



mode in the synchronous gauge) equal to zero is equiva- 
lent to simply taking the velocity component of the pres- 
sureless matter equal to zero (the comoving gauge of the 
pressureless matter); however, as will become clear in 
this work, even such a cure in the synchronous gauge 
is available only to the linear order in the presence of 
a zero-pressure fluid (Hwang & Noh 2006). The rem- 
nant gauge mode in the synchronous gauge causes quite 
a troublesome matter to handle in the nonlinear pertur- 
bation theory. 

Other gauge conditions free from the remnant gauge 
modes were proposed by Harrison (1967) for the 
zero-shear gauge (often known as the longitudinal or 
conformal- Newtonian gauge), and Nariai (1969) for the 
comoving gauge; see also Hawking (1966), Sachs & Wolfe 
(1967), and Field & Shepley (1968). There are, in fact, 
infinitely many other gauge conditions which are free 
from such a complication: i.e., free from the remnant 
gauge mode, and thus equivalently gauge invariant, see 
later. Systematic introduction of several different gauge 
conditions with explicit display of gauge-invariant com- 
bination of variables was made by Bardeen (1980) with 
huge success in later applications in the literature; see 
Kodama & Sasaki (1984) for a review. 

In a less known work, in 1988 Bardeen has suggested 
a pragmatic way of deploying the gauge conditions de- 
pending on the problems. As in other gauge theories the 
gauge choice is the degree of freedom which can be em- 
ployed depending on the advantages in achieving cither 
mathematical simplification or plausible physical inter- 
pretation. Bardeen has arranged equations so that the 
fundamental gauge conditions can be implemented easily. 
Bardeen's formulation of linear perturbation theory was 



2 



HWANG & NOH 



extended in Hwang (1991), Hwang & Noh (2001, 2005), 
and to the second order perturbations in Noh & Hwang 
(2004) and Hwang & Noh (2007). 

Our aim in this work is to extend the Bardeen's for- 
mulation to the exact and fully nonlinear order in per- 
turbations keeping the gauge advantages intact. We will 
display some applications, but since the main point is 
to present the new and powerful nonlinear perturbation 
equations, we will show detailed steps needed for the 
derivation in the Appendices A and B. The main equa- 
tions are presented in Section [3l and in the Appendix D 
for multi-component case. 

In Section [5] we review Bardeen's formulation and the 
gauge strategy in nonlinear perturbations. In Section [3] 
the exact and fully nonlinear equations are presented as- 
suming scalar- and vector-type perturbations of an ideal 
fluid in a flat background, but without fixing the tempo- 
ral gauge condition. In Section |4] we present equations 
valid to the third order in perturbations still without fix- 
ing the temporal gauge condition. In Section [5] we make 
several applications available in the comoving gauge in- 
cluding comparison with the Newtonian results. In Sec- 
tion [6] we consider the cases in other fundamental gauge 
conditions. In Section [7] we consider vorticity genera- 
tion from pure scalar-type perturbation in the comoving 
gauge; we will notice the gauge dependence of the con- 
cept of vorticity (rotation) to the nonlinear order. In 
Section |8] we present the equation of gravitational waves 
generated from pure scalar- and vector type perturba- 
tions which also depends on the gauge choice. In Sec- 
tion [3] we analyze the case of a minimally coupled scalar 
field. We show that in the comoving gauge the ideal fluid 
equations remain valid for the scalar field with a particu- 
lary simple equation of state. In Section [10] we comment 
on the possible future extension of this work. In the 
Appendix C we present analysis based on the covariant 
formulation. We set c = 1 . 

2. CONVENTION AND GAUGE STRATEGY 

Here are our metric and the energy-momentum tensor 
conventions. We consider scalar- and vector-type per- 
turbations in a flat Robertson- Walker background. The 
metric can be written as 

ds 2 = -a 2 (1 + 2a) drf - 2a 2 f/3 * + B$° A drjdx* 



+a 2 (l + 2ip)Sij 
xdx 1 dx :i , 



2 7,-, 



C 



(v) 



c 



2C. 



(t) 



(1) 



where a is the cosmic scale factor, and we assume 



i, and C [t) \ = = G {t)j 



*,3 



with 



dices of B\ , C\ v and Cy based on 5{j as the metric; 
indices (v) and (t) indicate the vector- and tensor-type 
perturbations, respectively. Indices a, b, . . . indicate the 
spacetime indices, and . . . indicate the spatial ones; 
we follow the convention of Hawking & Ellis (1973). 

The energy momentum tensor is given as (Ehlers 1993; 
Ellis 1971, 1973) 



T ab = tJ-U a U b + p (u a U b + g a b) + T^ab, 



(2) 



where \x and p are the energy density and the pres- 
sure, respectively, u a is a normalized fluid four-vector 



with u a u a = —1, and 7r a b is the anisotropic stress with 
7?ab = Kba, = 0, and TT a bU b = 0; tildes indicate covari- 
ant quantities. In this work we consider an ideal fluid 
with TT a b = 0. We set 

ju = fj, + Sfx = (i (1 + 6) , p = p + Sp, Ui = avi, (3) 

where /i and p are the background energy density and 
pressure, respectively. We set 



-Va + v 



(4) 



with v^ 1 i = 0; Vi and v\"' are based on 5tj as the metric. 
In this work, unless mentioned explicitly, we do not as- 
sume the perturbed metric and fluid quantities are small. 
The decomposition of an arbitrary spatial vector into 

(v) 

longitudinal and transverse parts as Bi = p^ + B\ , and 
a symmetric spatial tensor into longitudinal, trace, trans- 
verse, and tracefree-transverse parts as CV,- = ¥"%+7,ij + 

|(G|J + CjJ) + C[f are possible (York 1973); here, all 
spatial indices are based on Sij as the metric. We assign 
the transverse part as the vector-type perturbation, and 
the tracefree-transverse part as the tensor-type pertur- 
bation, and the remaining longitudinal and trace parts 
as the scalar-type perturbation. We can show that the 
decomposition is possible order by order in perturbation 
to nonlinear order. However, only to the linear order 
in the spatially homogeneous-isotropic background the 
three types of perturbations decouple from each other. 
To the nonlinear order we have couplings among the 
scalar-, vector- and tensor-types of perturbations. 

Here we ignore the tensor-type perturbation. Restric- 
tion of our attention only to the scalar- and vector-type 
perturbations can be regarded as our main assumption 
in this work. In Section [8l though, we will consider the 
contribution of tensor-type perturbation to the linear or- 
der. 

In considering the linear perturbation theory, Bardeen 
has suggested to take the spatial gauge condition 



,(«) 



7 = = CI 



(v) 



(5) 



but has saved the temporal gauge condition for later use. 
To the linear order 



X = a/3 + a 2 7, 



iC 



(«) 



(6) 



are spatially gauge- invariant, see below Equation (fig]); an 
overdot is a time derivative based on background cosmic 
time t (dt = adui). Under the spatial gauge condition in 
Equation (JSJ) our metric convention becomes 



ds 2 = -a 2 (1 + 2a) drj 2 - 2axidrjdx l 



+a 2 (1 + 2ip) 5 lj dx i dx\ 



where we set 



(8) 



We will take the metric and the energy-momentum ten- 
sor conventions in Equations ((TJ and ([2|) even in non- 
linear perturbation theory. Justification for taking the 
spatial gauge condition in Equation (JS|) to the nonlinear 
order will be made later in this section. It is essentially 
these spatial gauge conditions together with ignoring the 
tensor-type perturbation which allow our fully nonlinear 
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and exact formulation available. As will be explained, 
we do not lose any generality or convenience by taking 
these spatial gauge conditions; the only other alternative 
choice of spatial gauge condition leaves remnant gauge 
modes even from the linear order, see below Equation 
OH). 

The scalar-type perturbation equations without taking 
the temporal gauge condition are arranged in the follow- 
ing form by Bardeen (1988) 



ditions 





o • A 


(9) 


AttGS^i + 


A 

Hk+ -xcp = 0, 
a 1 


(10) 


A 

K + 


- 12nG(n + p)av = 0, 


(11) 


k + 2Hk 


- 4ttG (6(jl + 36p) + (^3H 4 


- 4) a = 0,(12) 


X + Hx- 


- ip — a = 0, 


(13) 



5fi + 3H (5(i + 8p) -(fi+p) ( k- 3Ha + -Av 



-0, 

[«VpR_i a 

a 4 (/Lt+p) a a([i+p) 



Sp 



(14) 
(15) 



where H = a/a; here we have assumed the flat back- 
ground and an ideal fluid. These equations are arranged 
without taking the temporal gauge condition. One ma- 
jor advantage of this arrangement is that the equations 
are designed so that we can readily impose various fun- 
damental gauge condition by simply setting one of the 
variables equal to zero. The vector-type perturbation 
equations without taking the temporal gauge condition 
are (Bardeen 1980) 



a 4 (fx + p) 



0, 



0. 



(16) 
(17) 



These vector- type equations are gauge- invariant. In this 
work we will present exact and fully nonlinear extension 
of these equations: see Equations (f2T|) -(|27 |) . 

The spatial gauge (congruence) condition in Equation 
(|5|) is a unique choice without remnant spatial gauge 
mode after taking the gauge condition. Thus, we do 
not lose any advantage upon our choice of the spatial 
gauge condition. On such a spatial gauge condition, the 
remaining variables can be regarded as spatially gauge- 
invariant ones (Bardeen 1988); see Equation ([6]). 

Concerning the temporal gauge condition, however, 
we have several fundamental gauge conditions available 
most of which remove the temporal gauge mode com- 
pletely. We have the following fundamental gauge con- 



comovmg gauge : 
zero — shear gauge : 
uniform — curvature gauge : 
uniform — expansion gauge : 
uniform — density gauge : 
synchronous gauge : 



V 


0, 


X = 


0, 


9 = 


0, 


K = 


0. 


S = 


0. 


a = 


0. 



(18) 



Also available as the gauge conditions are any non-gauge- 
invariant combination of these gauge conditions, thus we 
could have infinite number of different temporal gauge 
(spatial hyperspace or slicing) conditions available. As a 
consequence, we can manage arbitrary form of differen- 
tial equation for any non-gauge-invariant variable using 
a certain, perhaps ad hoc, choice of the gauge condition 
(Hwang et al. 2010). 

Bardeen's arrangement of equations apparently allows 
the simple adaptation of any fundamental gauge condi- 
tions: we simply set a perturbation variable equal to be 
zero. Except for the synchronous gauge the other tem- 
poral gauge conditions together with the spatial gauge 
condition in Equation (JS|) completely remove the gauge 
(coordinate transformation) degrees of freedom. Thus, 
each of the remaining perturbation variables has a unique 
counterpart of gauge-invariant combination involving the 
variable concerned and the variables used in the spa- 
tial and temporal gauge conditions; see below Equation 
(fH?|) . Therefore, all the variables in such a gauge condi- 
tion can be equivalently regarded as the corresponding 
gauge-invariant variables. 

The gauge conditions in Equations (0 and (|T8|) . and 
the above statements about gauge issue remain valid to 
fully nonlinear order as long as we take the nonlinear 
perturbation approach: this was shown in Section VI of 
Noh & Hwang (2004). In the following, we explain it 
again. 

We consider gauge transformation properties under 
x a = x a + C(x e ) with £° = £°, e = C, and & = 
£,,i/ a + Ci w ith ^ v ' i i = 0; C is based on S{j as the 
metric. To the linear order we have (Bardeen 1988; Noh 
& Hwang 2004) 



X = X 



v = v — £°, a = a , 

a y ' 

7 = 7- (p = f~ aH£°, 



B 



3H 



(v)l 



C. 



(tO 



a 
A 

TP 

(v) 



(v) 



(19) 



where a prime denotes the time derivative based on 77. 

Apparently 7 = = in all coordinates leaves = 0, 
thus fixing the spatial gauge degree of freedom com- 



(v) 

pletely; the only other choice taking p = = B\ ? in 

all coordinates gives £ ^ ^ , thus leaving remnant 
spatial gauge modes. Similarly, for the temporal gauge 
condition, for example, v — in all coordinates leaves 
£° = 0, thus fixing the temporal gauge degree of free- 
dom completely. The following combinations are gauge 
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invariant 

ip v = p — aHv = —aHv v , 

Vx = p-H x = -H X <p, (20) 

etc. This shows a systematic notation of expressing 
the various gauge-invariant combinations. This nota- 
tion is practically useful to implement the spirit of 
Bardeen's formulation employing many gauge conditions 
which make all variables gauge invariant (Hwang 1991). 
The gauge-invariant combination, for example, tp v is the 
same as ip in the v = hypersurface condition, thus 
ip v = ip\ v =Q. The temporal gauge condition, for exam- 
ple, v = fixes the temporal gauge mode completely. 
Thus, any perturbation variables in that gauge, for ex- 
ample ip, can be equivalently regarded as a temporally 
gauge invariant ones, i.e., ip\ v — <p v . Similar complete 
gauge fixings are true for the other temporal gauge con- 
ditions. The synchronous gauge (a = 0) is an exception, 
leaving a remnant temporal gauge mode £ (?7, x) oc a -1 
even after fixing the gauge condition. 
Now, to the nonlinear order, we may set £° = f ' 1 ) + 

£°( 2 ) + . . . and & = Q 1 ' + ^ + . . . where the number 
inside the parenthesis indicates the perturbation order. 
To the second order, the gauge transformation properties 
of each variable have the same form as in Equations (|19[) 
with additional terms involving quadratic combinations 
off , £ and perturbation variables, all to the linear order. 
Since each of the quadratic terms involves £ or £ to the 
linear order, as long as we take the spatial and temporal 

gauge conditions which lead to = = fj 1 ' (thus 

the synchronous gauge is excluded) , we have exactly the 
same gauge transformation properties in Equation (| 19[) 



now valid for pure second order variables: for example, 
we have — — S, ^, etc. Thus, by imposing the 
same (i.e., ones removing the gauge degrees of freedom 
completely) gauge conditions now to the second order, 

we can remove f°( 2 ) = = thus leaving any vari- 
able in that gauge having a corresponding unique gauge- 
invariant counterpart: i.e., ip\ v — (p v , etc. Apparently, 
the same process can be continued to any higher order 
perturbations. Explicit forms of the gauge transforma- 
tion properties and gauge-invariant combinations to the 
second-order perturbation are presented in Noh & Hwang 
(2004) and Hwang, et al (2012). 

The names of our gauge conditions can be justified 
to the nonlinear order by examining the ADM met- 
ric, extrinsic-curvature, intrinsic-curvature, and the fluid 
variables presented in the Appendices B and C. For the 
comoving gauge with v = 0, ignoring the vector- type 
perturbation, from Equation (|110|) we have Ui — 0, thus 
the fluid four-vector becomes the normal four-vector. For 
the zero-shear gauge with \ = 0, ignoring the vector- type 

perturbation, from Equation (|109[) we have if* = 0, thus 
having vanishing shear of the normal flow vectors n a ; we 
have 3y = —Kij, see the Appendix C. For the uniform- 
curvature gauge with ip = 0, from Equation (|108[) we 
have R jj^ = 0, thus having vanishing curvature of the 
spatial hypersurface; in the presence of the background 
spatial curvature, we have spatially uniform curvature. 
For the uniform-expansion gauge with k = 0, from Equa- 
tion (|109p we have K = —3H, thus having the trace of 
extrinsic curvature uniform; we have 8^ =n c . c = —K, 
see the Appendix C. For the uniform-density gauge with 
5 = 0, from Equation ([3]) we have Jl — fx, thus the density 
becomes uniform in the hypersurface. 



3. EXACT AND FULLY NONLINEAR PERTURBATION EQUATIONS WITHOUT TAKING TEMPORAL GAUGE CONDITION 

Fully nonlinear extension of Equations (|9"j)- (fTT)) will be presented below. These equations are the main result of this 
work. Based on the ADM equations, the derivation of our fundamental equations is unexpectedly simple. In order 
to help the reader who will attempt the derivation we review the ADM formulation in the Appendix A, and present 
detailed steps required for the derivation in the Appendix B. 
Definition of n: 



1 



1 



Mil + 2<p) 



3ip - 



1 



x k , k 



l + 2<p 



(21) 



ADM energy constraint: 



BttG. 



A 

—^3 



Hk 



-tt, — — tt: = — K 2 — 4:7rG (u + p) 1 'Z 7 t Pj 

a 2 (l + 2p) 2 6 y ^ Fl l + 2ip 2a 2 {l + 2pf 



(22) 



ADM momentum constraint: 



1 



2aW(l + 2<p) 
1 f 



AXi + 3 X\ ik ) + (/* + p) aviJ 1 + - 



V k Vk 



2tp 



Ma 



a 2 M(l + 2p)\ \ M 



l + 2pj [2 V 



xi 



<p>- 



(i + 2py 



M „ 

l + 2<^ 3 



"77 ( X°f,i + Xi 1 ? 3 ~ -Sfx k f,k 



(23) 
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Trace of ADM propagation: 



1 



3— -.H - 3H Z - 4ttG (p + 3p) + A + — -.k + 2Hk 



1 



a 2 7V(l + 2ip) 



x « i + 



l + 2tp 



Tracefree ADM propagation: 



1 d 



X 



1 



a 2 N{l + 2cp) 
1 



^(x\,+x/)-^A- T ^-(xV J + 



a 2 {l + 2<p) 
V- + V 



1 



1 + 2^3 



3 J 



A/" 



= 8nG 



l + 2<p 



1 



2^ 
1 



1 



Xj<P'' 



aW 2 (l + 2(^) 2 



1 2 

+ 1 + 2 ^ {x' l XkV,j - X h XjV,k + XfcjxV' - Xi.fcXV ) + (x + 2y ,)2 b^XiV <P,k - X XfcV'V.j) 



1 



02(1 + 295)2 

ADM energy conservation: 

1 /a 



1 + 2^ 



r 



N \dt a 2 (l + 2</j) 



M + + p) 



1 + 2<^ jvy 0(1 + 2^) 



1+2^ 

/x + p 



+ (j5 + p) (3ii - k) 1 + 



4 irv fe 
3 1 + 2(f 



v\ 1 + 



V k Vk 

l + 2<f 



(24) 



(25) 



M + P 



a 2 AA(l + 2(^)2 
ADM momentum conservation: 



1 2 

Y ,J VjVi Y „-l> 1>,- 

A 3 3 A J l + 2ip 



X J 



1 d 



a (fj,+p)vi\ 1 + 



V k Vk 

l + 2<p 



l + 2tp 



(26) 



(^ + p)u 3 -Wl + 



v k Vk H+P j 

l + 2tp ~ 1 + 2^ U 



(27) 



Equation (|2 1|) follows from the definition of « as K = —3H + k; K is the trace of extrinsic curvature presented in 
Equation (fT09]l , Equations p2]) -([27 )l follow from the ADM equations in Equations ([5? ]) - (jT02"]) in the Appendix A. The 

variable M is related to the lapse function in Equation (|105l) . and if* is the tracefree part of extrinsic curvature in 

Equation (fTU5j) . With N and kM\ given as 



N= a/1 + 2c 



X fe Xfe 



a 2 (l + 2^)' 



1 



1 + 2^ 



\x>^ (x»,j + Xi,i) - \x l ti X 3 f>,3 



aW 2 (l + 2^) 2 
2 



-jX l,J (Xi,j+Xj,i) - \x\ix\ 



1 



(1 + 2^)2 



2 X XW J< P,i + oX XV,i^,j 



(28) 



Equations (|2ip ~ (|27p are complete set of fully nonlinear perturbation equations valid for the scalar- and vector-type 
perturbations assuming ideal fluid in a flat background; instead of Equations (|26p and (|27[) we can use alternative forms 
in Equations (|124l) and (|125[) . respectively. We can consider this set of equations as exact, or treat it perturbatively 
to fully nonlinear order. The perturbation variables are 5fi 7 Vi, K, Xi> f an d ct; Sp should be provided by an equation 
of state. Notice that we have not separated the background order equations yet; we only have assumed that a is a 
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function of time. The vector-type perturbation is contained in Uj and Xi as 

Vi = -v ti + vl v) , Xl = X,i + a^[ v) ■ (29) 

For the vector- type perturbation, Equations ([2"B")) and (|2T|) to the linear order give Equations (IT51) and (jTTJ) , respectively. 

For the pure scalar-type perturbation we set v\ v ^ = = ^ , thus V{ = —vj and %i = x,i- 

In the above set of equations we have not taken the temporal gauge condition yet. In a sense the equations are 
in a sort of gauge- ready form. As the temporal gauge condition we can impose any one condition in Equation (|T8)) , 
except for the synchronous gauge which leaves the remnant gauge mode; see explanation in the next paragraph below 
Equation ([20]). Thus, as the gauge conditions we have 



comoving gauge : v = 0, 

zero — shear gauge : X = 0, 

uniform — curvature gauge : ip = 0, 

uniform — expansion gauge : k = 0, 

uniform — density gauge : 5 = 0, (30) 



now valid to all perturbation orders. These are the fundamental gauge conditions available to fully nonlinear order. 
Also available as the gauge condition is setting any non-gauge-invariant linear combination of these fundamental gauge 
conditions equals to zero. We can also take the different gauge condition for the different perturbation order. In these 
ways we have infinite number of gauge conditions available, which was true even to the linear order, now to each 
perturbation order. Under these gauge conditions which remove the gauge mode completely, all perturbation variables 
have unique gauge-invariant counterparts, thus we can identify these as the gauge-invariant variables. Therefore, 
the nonlinear perturbation variables in any of our suggested gauge conditions mentioned above can be regarded as 
gauge-invariant ones. 

In the Appendix D we will present the case of multiple-component fluid system. The above set of equations remain 
valid with the fluid quantities interpreted as the collective ones. We will present the relations of the collective fluid 
quantities with the individual ones, see Equations (|132[) and (|137j) . and present the additional energy and the momen- 
tum conservation equations of the individual component, see Equations (|143j) - (|146j) and the prescription explained 
above Equation p43[) . 



4. THIRD-ORDER PERTURBATION EQUATIONS IN A GAUGE-READY FORM 

In the nonlinear perturbation approach we assume the perturbation variables 5, Vi, k, Xi, ot and cp are small. As one 
example, here we present pure scalar-type perturbation equations valid to the third order in perturbations without 
fixing the temporal gauge condition, thus Vi — —v t i and Xi — X.i- Equations up to third-order perturbations are 
needed to get the leading nonlinear contribution to the power spectrum. 
Definition of k: 



3Ha 
3 



A 



3ip 

Ha 2 (-3 + 5a) 



X X,k [3H + 3Lp - 6Hip - 9Ha - 
A \ f _ 3 



A 



rX'V.fe (1 - ot - Acp) 



'Sip 



rX 



2tp a - 2atp - Atp' 



(31) 



ADM energy constraint: 



■I ' " J 



8ttG A\ A 
-/i — — ] + inGSfi + Hk + —ip 



-k 2 - 4irG (fj, + p) v>*v i I 1 - 2cp ■ 
6 



a 

<5/i + Sp 
M + P 



I ( ^<fi) <P (1 - ty) + (1 - M 



1 

~4^4 



x l3 x,a 



i(Ax) ! 



(1 - 2a - 4ip) - ±X' i3 X,iV,j + o ( A X) rV,< 



(32) 
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ADM momentum constraint: 



7 



k + —x - 12ttG (fj, + p) av 



a 

2 / A 

3 W 



= 8nG (ji + p) av.i ( ^J^y + \v\k) 



3 1 

a + 2ip - -a 2 - 2a V - 4ip 2 + 7^2 X' fc X,fc 



1 



a ,j (1 - 3a - 2<p) - ip j (1 - a - 4p) + -^X' k X,kj 



+ - 



(A X ) <p,i + x,;A<^ + - {r k V,k) ti 



(1 - a - Aip) 



XW k (" fc + <p,k) + gX ,fe (<P,kV,i + 3^,,a,fe - 2<p )fe a,j) 



(33) 



Trace of ADM propagation: 



3H - 'AH 2 - AttG {h + 3p) + A + k + 2Hk - 4irG [dp + 35p) + [ 3H + ) a 



■ ( - ■ - + ^X' k X, k ) +\k 2 + 8nG (p + p) v\i (l-2 V + ^±^) 



1 



3a 2 - ^X k X,k (1 - 3a - 2 V ) - 5a 3 



a + 2cp- \ 2 - 2aip - V + ^X' fe X,fc J 



(1 - a - 2ip) 



A 

2^2 



a 2 - \r k X,k (1 - a - 2<p) - a 3 



" ^2 [x' </s .< (1 - a - 2p) + p-'a.j (1 - 2a - Asp)] 



X ' V XM - 3 (Ax)^ 



(1 - 2a - 4y>) - 5x ,y X,»V,j + 3 (A*) X'V, 



(34) 



Tracefree ADM propagation: 



- ( V V, --6*A\ (x + H X -a-<p)=ll + Ml)\ - 



1 

a 2 



a + 2ip--a 2 - 2aip - 4tp 2 + — X ' k x, k 



2a 2 ' 



2 „, 



x ( V 4 V, - |^A j x + I x'V J + XjV' 1 - §<Jjr V* ) (1 - a - 4y>) 



(1 - a - 2<p) ( V'V, - -«5;-A ) x 



2 



-1 (a + 2^ - |a 2 - 2a^ - 4^ 2 + ^X'^X.fc) (Vv, - 1^ a 



-a 2 + ^r k X,k (1 - a - 2ip) + a 3 



(v'Vj - ^8j<P' k <P,k) (1 - M - ^2 (''",; + VJ"'* - (1 - 2a - 4<^) 

~<p (1 - 3p) (Vv, - ^}a) p + SttG (^ + p) (V^- - 1( 



<5/z + 
p + p 



(35) 
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ADM energy conservation: 

fi + 3H (fx+p) + 5(x + 3H (5(i + Sp) — (fj, + p) k — jxa — (fx + p) —v = (Six + Sp) k + ^ (ix + p) nv' k v k 

a " 3 



1 



i 



-3a 2 + -Vr fe X,fc (1 - 3a - 2tp) + 5a 3 



1 



1 



Sp [ a - -a 2 + X,k ) - ^X'^M.i (I -a- 2tp) 



+-(p + p) V^v 
a 1 

+- (p+p)v' i ( 1 + 



V M + P / 

<5/i + Sp 



-4JI (p,+p)v le v, k \ l-2(p 



5/i + Sp 
M + P 



-2ip- 

Sji + Sp 
fi + p 



H + p 



(l-2p) + 4<p 2 + -v> k v, k 



[3ip,i (1 - 4ip) + 2a,, (1 - 2a - 2<p)] 



+ —{p + p)v l 



2 1 



(36) 



ADM momentum conservation: 
1 



[a 4 (p + p)v\ — (fx+p)a — Sp 



Of 



3H 



a(fx+p) v., 



Six + Sp 1 k 
M + P 2 



3 2 1 fc 

a_ 2 a + 2^ X ' X ' fc 



dt 

Sp + Sp 
P + P 



:i-a-2(p)x' fe V fc 



a (/i + p) v.i 1 + 



Sp + Sp 
P + P 



(fi+p) V,k I 1 + 



5/i + (5p 
M + P 



X' K 4 (1 - a - 2y>) - 2 X > h <p 



+ ([i + p)v k v }i (a + 3ip) k +v j k(a-ip) 
+ (fi+p)\a 



2a+ S_E±Sp (1 _ 2a) + w 



M + P 



1 

2^2 



(5/i + Sp 
P + P 



\ [x' fc X,fe(l-2a-2 ¥3 )] ! 



To the background order, Equations (f3"2"j) . (jM)) and (|36p . respectively, give 



H z = 



8ttG A 



3 3 a 



^(/i + 3p) + ^, ft + 3H(p + p) = 0. 



(37) 



(38) 



To the linear order, assuming the background equations separately, Equations (j3"Tj) - (|37|) give Equations ([9|)- (fT5|) . 



5. COMOVING GAUGE 



In sections [5] and [6] we consider only the scalar-type perturbation, thus Vi — —vj and Xi = X,i- I n the comoving 
gauge we set 



v = 0. 

We have Uj = 0, thus the fluid four- vector becomes a normal one with Ui = 0. Equation (f2~T|) gives 

P,i = - (m + p) 



(39) 



(40) 



Equations (j2"l~j) - ([2"7|) give a (redundantly) complete set of equations for the variables S, k, <p, x an d a. We can treat 
this set of equation either exactly or perturbatively to all orders. As explained in Section [5] and below Equation (|3T)]) . 
all perturbation variables in the comoving gauge are gauge invariant to the nonlinear order as 



S = 5 V , k = k v , (f = ip v , X = Xv, a = a v . 



(41) 



We have several ways of having closed form second order differential equations. (I) Equations (f2l)f and (|M|) . together 
with Equations (f2"2"j) . (|2U1) and (|27p to determine a, y and x, give equations for S and /i. (II) Equations ([2"T]) and (f2"4")l . 
together with Equations (f2"2"l) . ([2"B"f and ([2"7]) to determine <5, a and x, give equations for <p and «. Etc. For example, 
assuming the background equations are valid separately, Equation (|26p gives 



1 1 (d 



p+pM\dt a 2 (l + 2ip) 



Vi U 



/i 



M + P 



(42) 



FULLY NONLINEAR COSMOLOGICAL PERTURBATIONS 
By removing k in Equation (|24p we have a second-order differential equation for S as 



1 

77 



i l 



d 



X 



Jl + pAf\dt a 2 (l + 2if) 
1 

Jf 



/' 



fi+p 



2H 



1 1 



d 



X 



Jl + pJ\f\dt a 2 {l + 2ip) 



V, fJL 



H + p 



1 if - 4ttG (6/j, + 3Sp) 



AJ\f 



aW(l + 2^) 



aW(l + 2y>) 



X «,< 



l + 2<p 



(43) 



Similar equation can be derived from Equation (|123j) evaluated for Uj = 0. 

If we have a solution for a given variable we can derive the rest of the variables using the above complete set of 
equations. From these solutions we can derive all the variables in any other gauge conditions using the gauge-ready 
form equations in (|21[) -(j27 p and explicit construction of gauge- invariant combinations; for solutions in the matter 
dominated era to the second order perturbations, see Hwang, et al (2012). 

5.1. Zero-pressure case 
For p = 0, Equation ([4TJ]) gives A/"j = 0. Thus, we can set Af = 1 with 



1 X k X,k 



2a 2 (l + 2f)' 



(44) 



Thus, in our spatial gauge condition with 7 = 0, the comoving gauge (v = 0) is no longer synchronous (a = 0) to the 
nonlinear order even in the zero-pressure medium (Hwang & Noh 2006). 

Equations and (|26l) together with Equations (|2"2"|) and (|2"3")l provide a complete set of equations for S and k. We 
have 



-+3H)(1 + S)+S-k = Sk- 



a 2 {l + 2ip)' 



1 



-3H - 3H 2 - AttG^l + A + k + 2Hn - 4ttG<5/x = -k* - 



X «,f 



a 2 (l + 2<p) 



(45) 
(46) 



with % and y determined by 
Ax,* 1 



a 2 (l + 2(p) a 2 (l+2ip) 2 



2(A X )^: 



2 

where 



BttG 



A 

—^3 



Af 



a 2 (l + 2ipf 



3 3 1 

2 X ' 1 ^ ~ 2 1 + 2<^ 



<P' <P,i 



2a 2 (l + 2tpY 



1— i— j 



(47) 
(48) 



a 4 (l + 2<p) 5 



x < X' J X,-. 



1,4 



3 (Ax)x'V.t -X %3 X,if,j 



(1 + 2^ 



1 2 

3 (x'V,i) +x ,, x, l vV,; 



(49) 



These equations are still exact. Assuming the background equations are valid separately, from Equations (|45[) and 
(|46|) together with Equations (j47|) and Equation (|48l) . we can derive a closed form second-order differential equation 
for 8 or n. For <5, we have 



6 + 2H6 ~ AnGfiS 



1 



(a 2 KS)' 



1 



X% 
1 + 2^3 



1 2 r% 
3 K ^ 



a 2 (l + 2v?) 



f^. 



(50) 



Using Equations (gSJ, gTJ) and 
terms of <5. 

Now, equation for ip follows from Equations (f2"Tj) and (|4T|) as 



to determine k, tp and Xj we can express (perturbatively) this equation purely in 



[In (l + 2p)V 



1 



a 2 (l + 2^) 2 



X' V.ife + X,i Al P 



1 



1 + 2(^9 



(51) 



thus ip — to the linear order. Together with Equation (|46[) . using Equations (|47|) and (|48| to determine 5 and x, we 
have the closed form equations purely in terms of ip and k. 
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5.2. Equations to fifth order in the zero-pressure case 

As one exercise demonstrating the power of our fully nonlinear formulation, we present the perturbation equations 
to the fifth order. We consider a zero-pressure fluid in the comoving gauge. Obviously the derivation is quite simple 
requiring only Taylor expansion of 1/(1 + 2ip) terms. Assuming the background equations are valid separately, from 
Equations (|35 ]) -(|jg ]) we have 



k = Sk - -rY'** i (l - 2w + 4ip 2 - 8ip 3 



1 



1 



k + 2Hk - AnGSfi = -k o X'X* (1 - 2<p + 4tp z - 8<p" 



k)kI 



(52) 
(53) 



2 ( 

— K 

3 I 


+ H, 


1 

~ 2 


H 


- [xWa - 


I" X,i<P' 3 


Hk + 


AirGfiS + 


A 

— V = 



- (A X ,i) V (1 - + V - 8^ 3 ) - X ,3 i - ~<^A X <Pj (1 - 4^ + 12^ 2 - 32^ 



+ (x' j V,i + X,if' J ' - \4x' Vfc) (1 - 4^ + 12^ 2 - 32<^ 3 ) - (x,if,j + J\,r-.) </> J (l-6<p + 2 V) 



(54) 



A 



3 1 
2a 



-^V i (i-6p+2V 



3 ) - ^X(55) 



where 



+2 



1 2 

3 (x'V.i) + X'^iVVj 



(l - Aip+ 12tp 2 - 32v? 3 ) 
(l-8p) L 



(Ax) X'V.i - X ,ij X,i<P,; 



(l-6<p + 2A(p 2 



(56) 



From Equations (1521 and (JS3J) together with Equation (1541 to the fourth order and Equation ([5"5j) to the third order, 
we can derive a closed form second-order differential equation for S or k. 

The fifth order perturbation equation is needed to have the next-to-leading-order nonlinear contribution to the 
power spectrum. The leading nonlinear order power spectrum demands third order perturbation and the results for 
the density and velocity power spectra are presented in Jeong et al (2011). The fourth order perturbation equations 
will be needed to have the next-to-leading-order nonlinear contribution to the non-Gaussianity. The leading order 
non-Gussianity demands second order perturbation. 



5.3. Pure Einstein's gravity corrections to fully nonlinear order 
We can arrange Equations (|4"S"|) , (|4*6")) and (|4"7) in the following forms 



5 — k — 5k H — ^X l °~ i 

k + 2Hk - AttGSh - 
2 

+ 



a 2 {l + 2<pY 

1 2 1 ^ 1 

K + — X K>i - 



r ij x,a - 1 (Ax) 2 



2<px' t K. 



Acp{l + if) 



a 2 (l + 2ip) a i (l + 2tp) 2 



a 4 (l + 2(p) 3 
A \ 2<pA x , t 

2* j „2 



-(A x)r V )i -2x-X,^ + TT ^ 



1 



(x'V,i) 2 + X' l X,i ( P' 3( P,: 



a 



a 2 (l + 2^) 



+ 



1 



a 2 (l + 2^) 2 



1 



1 3 

2 (Ax) tp,i + ' <P,ih ~ X,ik<P' + -jXAf - 2 , _ 



(57) 



x' l3 x,a 



i(ix) 2 



X,i<P,k + ^X,k<P,i ) V 



(58) 



(59) 



Terms in the right-hand-sides are pure Einstein's gravity corrections; see the next Section. In the zero-pressure 
case Newtonian perturbation is closed at the second order in perturbations (Peebles 1980; Vishniac 1983; Zel'dovich 
& Novikov 1983). Notice that the pure Einstein's gravity contributions involve <p which is related to the (spatial) 
curvature perturbation in the comoving gauge, see Equation (llOSp . 



5.4. Relativistic/ Newtonian correspondence 

Except for terms involving tp Equations (|57 j) -(|59 j) co- 
incide exactly with the Newtonian hydrodynamic equa- 
tions of the mass and the momentum conservation equa- 



tions (removing the gravitational potential in the mo- 
mentum conservation equation using the Poisson's equa- 
tion), respectively. This statement is true to fully non- 
linear order in perturbation in the presence of the cosmo- 
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logical constant in the background. That is, by ignoring 
p terms (we cannot do this in general, though) Equations 
gU), gU) and gZl) give 



5k 



2Hk 
1 



1 ,ix 

2* °- 

or 

4ttG (Sfj, - 
1 , 



36p) 



K.: 



a- 



■: ij x,a - \ (Ax) 2 



A 

— x = o. 



(60) 



,(61) 
(62) 



By identifying 8 and u as the Newtonian density and 
velocity perturbations with 



1. 



thus x — au with u = Vu, we have 



(63) 



5+-S7 ■ 
a 

1 

-V- fu- 
el 



1. 



u = — V • (Su) 

a 



(64) 

Hu) + AnGgS = — • (u • Vu) , (65) 



1 



which coincide exactly with the Newtonian hydrody- 
namic equations with the Newtonian gravitational po- 
tential removed (Peebles 1980; Vishniac 1983; Noh & 
Hwang 2004). 

Therefore, to the fully nonlinear order the pure Ein- 
stein's gravity contributions appear in terms of ip. The 
ip terms start appearing from the third order pertur- 
bation, thus we have exact relativistic/Newtonian cor- 
respondences of the density and velocity perturbations 
in the comoving gauge to the second-order perturbation 
(Noh & Hwang 2004). In the comoving gauge, however 
we do not have correspondence with the Newtonian grav- 
itational potential. In the conventional Newtonian limit 
of Einstein's gravity, a usually corresponds to Newtonian 
gravitational potential. As in Equation (|44p . a vanishes 
to the linear order, and we no longer have the proper 
Newtonian correspondence for the gravitational poten- 
tial in the comoving gauge. 

6. OTHER FUNDAMENTAL GAUGES 

6.1. Zero-shear gauge 

In the zero-shear gauge we simply set 

X = 0. (66) 

We have iV* = = . The metric becomes 

als 2 = -a 2 (1 + 2a) drf + a 2 (1 + 2cp) 5 ij dx i dx j . (67) 

In linear theory, the zero-shear gauge is quite popular in 
the literature (Mukhanov et al. 1992), despite its short- 
comings in numerical treatment in the early universe and 
in the Boltzmann code (Ma & Bertschinger 1995; Hwang 
& Noh 2001). Equations (|2T ]) -(|27 1) give a (redundantly) 
complete set of equations for the nonlinear perturbation 
variables 5, v, k, p and a. All perturbation variables in 
the zero-shear gauge are gauge invariant with 

S = S X , v = v x , k = k x , <p = p x , a = a x .(68) 

There are several ways of having the closed form second 
order differential equations. (I) Equations (|26|) and (|24|) . 



together with Equations (|2"21 , ([2"B"f and (|25l) to determine 

a, ip and Xi gi ye equations for 8 and k. (II) Equations 
P6l) and (|27)l, together with Equations (US]), and 

(|23[) to determine a, p and k, give equations for 8 and 
v. Etc. 

It is known that in the small-scale limit (inside the 
visual horizon) in the matter dominated era, the den- 
sity, velocity and gravitational potential variables (8, v 
and a) in the zero-shear gauge coincide exactly with the 
Newtonian results (Hwang & Noh 1999); this is true 
even for the variables in the uniform-expansion gauge. 
We have shown that the correspondences continue to be 
valid even to the second-order perturbation (Hwang, et 
al. 2012). Whether this relativistic/Newtonian corre- 
spondences would continue to higher order in perturba- 
tion or valid in to fully nonlinear order is an interesting 
subject for future investigation. 

6.2. Uniform-curvature gauge 
In the uniform-curvature gauge we set 

ip = 0. (69) 



We have R {h)i jM = 
The metric becomes 



0, thus flat in the flat background. 



ds 2 = -a 2 (1 + 2a) dr] 2 - 2ax,idr / dx i + a 2 8 ij dx i dx j .(70) 

Equations (|2"Tj) - (l2~T|) give a complete set of equations for 
variables 8, v, k, x an( l a - All perturbation variables in 
the uniform-curvature gauge are gauge invariant with 

8 = 8 V , v — v v , k = k v , X = Xf: a = a v .(7\) 

In linear theory, the uniform-curvature gauge is use- 
ful to handle scalar field perturbation (Field & Shep- 
ley 1968; Lukash 1980a, 1980b; Sasaki 1986; Mukhanov 
1988; Hwang 1994; Hwang & Noh 2005). To the linear 

order, with the scalar field 4> decomposed as cf> = cf> + Scf>, 
we have S(f> — S(f> — </>'£°, thus 



H' 



(72) 



is gauge invariant. For fully nonlinear treatment of the 
scalar field perturbation in the uniform-field gauge (8<p = 
to the nonlinear order), see Section |H1 



6.3. Uniform- expansion gauge 
In the uniform-expansion gauge we set 

K EE 0. 



(73) 



We have 6^ = —K = 3H, thus the expansion rate of the 
normal frame vector field is uniform in space. Equations 
(|2"l"j) -(f2"7 | give a complete set of equations for variables 8, 
v, x, if and ct. All perturbation variables in the uniform- 
expansion gauge are gauge invariant with 



S = S K 



X = Xk, a = a K .(74) 



As in the zero-shear gauge, the uniform-curvature gauge 
also shows small-scale Newtonian correspondence of 8, v 
and a up to the second order in perturbation (Hwang et 
al. 2012). 



12 



HWANG & NOH 



6.4. Uniform- density gauge 
In the uniform-density gauge we set 
5 = 0, 



(75) 



thus the density is uniform in the hypersurface. Equa- 
tions (|2"Tj) - ([2"T| give a complete set of equations for vari- 
ables v, k, x, f and a. All perturbation variables in the 
uniform-density gauge are gauge invariant with 



vs, 



Kg, 



<P = V5, X = X8, a = a s . (76) 

Curvature perturbation <p in the comoving gauge ((p v ), 
in the uniform expansion gauge {(p K ) and in the uniform- 
density guage (ips) are known to be conserved in the 
large scale (super-sound-horizon scale) to the second or- 
der in perturbations (Hwang & Noh 2007), and to general 
nonlinear order based on the spatial gradient expansion 
method (Lyth et al. 2005). The proof based on our ex- 
act and fully nonlinear perturbation equations deserves 
a close examination. 



7. 



GENERATION OF VORTICITY IN THE COMOVING 
GAUGE 



We consider the generation of vorticity (rotation or 
vector-type perturbation) from the pure scalar-type per- 
turbation. We consider the comoving gauge. Thus, we 
have Vi = in the v% terms multiplied by perturbation 

(v) 

terms, and have Vi = v\ in the pure v% term without 
perturbations multiplied. From Equation (|27|l we have 

-P,-(P + P)^. (77) 



1 



a (fi + p)v 



Thus, we have 
1 



v) 



(78) 



which is true to the fully nonlinear order. 

In terms of the covariant equations, from Equation (8) 
in Hawking (1966) we have 



~r d 
i b 



Wed 



\c;d] 



= --8u>ab ■ 



2? c [Q S b]c , (79) 



where h a b = g a b + u a Ub is the projection tensor; for the 
covariant notations see the Appendix C. For the vorticity 
generation from pure scalar-type perturbation in the co- 
moving gauge, from Equation (|117[) we have Wij = av 



0) 



and from Equations (f79|) . (|12ip and ()122|) we arrive at 
Equation ([75| . As commented below Equation (|119p . 
to the nonlinear order the vector-type perturbation is 
no longer gauge-invariant and the above conclusions are 
valid only in the comoving gauge. This issue deserves 
further investigation. 

8. GENERATION OF GRAVITATIONAL WAVES FROM 
SCALAR- AND VECTOR-TYPE PERTURBATIONS 

Fully nonlinear and exact formulation including the 
tensor-type perturbation is supposed to be a complicated 
subject which is left for future investigation. Here we 
consider a much simpler case with linear tensor-type per- 
turbation. In the presence of the tensor-type perturba- 
tion a change occurs in the spatial part of the metric in 
Equation (|103[) as 



9ij = ° 2 [(1 + 2 V) $ij + 



(80) 



where hij is the transverse (/i^. = 0) and tracefree 

(h\ = 0) tensor-type perturbation; indices of hij are 
based on 7^ as the metric; only in this section hij in- 
dicates the tensor-type perturbation. By keeping only 
linear order terms in hij we can update quantities in the 
Appendix B. In our basic perturbation equations in (|21[) - 
(|27[) . the linear tensor- type perturbation contributes only 
in Equation (|25p by simply adding the following term 



n2 ' L V 



(81) 



in the left-hand-side. We ignore the tensor-type 
anisotropic stress. Let us write Equation (1251) including 
the linear tensor-type perturbation as 



A 

3Hhij ^ hij — n 



> j ■■ 



(82) 



where we moved all the terms in Equation (1231) to the 
right-hand-side and called it n*, then lowered the in- 
dex. The right-hand-side of this equation includes linear 
parts of the scalar- and vector-type perturbations which 
need to be removed to get the pure tensor-type pertur- 
bation equation generated by the nonlinear scalar- and 
vector-type perturbations. By the following operation 
we can separate the linear part of the scalar- and vector- 
type contributions [see Equation (210) in Noh & Hwang 
(2004)] 



. A 

hij -\- SHhij - 2^*J — ^ij) 

Sij = Uij - 2A~ 1 V( i n* : ))fe 



A- 2 (V^ + ayAJnV 



(83) 
(84) 



This can be regarded as the equation for gravitational 
waves (tensor-type perturbation) generated from pure 
scalar- and vector-type perturbations to the fully nonlin- 
ear order. The nonlinear terms in Sij still depend on the 
temporal gauge condition, and consequently the gravita- 
tional waves generated from the scalar- and vector-type 
perturbations do depend on the temporal gauge choice. 
In the perturbation approach hij should be regarded as 
the same order perturbation as the one considered in ; 

i.e., with an expansion hij = h^ + h^ + . . . , where in- 
dex (1) and (2) indicating the order of perturbation, for 

(2) 

quadratic order perturbations h^ is the same as h\j , 
etc., thus depending on the gauge choice. This subject 
deserves further studies. 



9. SCALAR FIELD IN THE COMOVING GAUGE 
We consider a minimally coupled scalar field <f> with 



9ab, 



and the equation of motion 



(85) 



(86) 
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Compared with Equation the fluid quantities can be 
read as 

p = f ab u a u\ p=^f ab h a \ 

nab = T cd h c a h d b -phab- (87) 

The fluid quantities and the equation of motion to the 
fully nonlinear order can be derived in a gauge-ready 
form. Here, we only consider a fluid formulation of the 
scalar field in the comoving gauge. From Equations ([2]) 
and ([55)) we have 

Ui = -lf lb u b = (88) 

H (j) yC U c 

Thus, the comoving gauge (v = 0) implies the uniform- 
field gauge {5<f> = 0) and vice versa to the fully nonlinear 

order; we set 4> = </ ) + ^<^ where <f> is the background order 
scalar field. In this gauge we can show 

5 = 2A^ + V ' P = W*^ ~ Vl * ab = °' (89) 
Thus, to the background order, we have 

V=\tf + V, p^^-V, (90) 
and to the fully nonlinear order, we have 

5p = s » = -ir^{ 2a+ ^Skr))> (91) 

with vanishing anisotropic stress. Therefore, the ideal 
fluid equations in Equations (|2T1) - (|27)) under the comov- 
ing gauge remain valid with the perturbed equation of 
state given as 5p — <5/i. 



10. DISCUSSION 

Extension or feasibility of similar fully nonlinear for- 
mulation including the following cases deserves future 
investigations: (i) background spatial curvature, (ii) 
anisotropic stress, (iii) a minimally coupled scalar field, 
(iv) multiple components of fluids and fields, (v) class of 
generalized gravity theories, (vi) the electric and mag- 
netic fields, (vii) the covariant equations and the Weyl 
tensors, (viii) the tensor-type perturbation, (ix) null 
geodesic for Sachs- Wolfe effect, (x) gravitational lensing, 
(xi) Boltzmann equations for photons, and massless and 
massive neutrinos, (xii) gauge transformation properties, 
(xiii) expression of gauge-invariant combinations, (xiv) 
equations in mixed gauge conditions, etc. At the mo- 
ment the first seven are trivial (some could be tedious 
though) extensions while the remaining ones need closer 
examinations for their feasibilities. Implementations of 
all the above cases (except for x) were made to the second 
order in perturbations in Noh & Hwang (2004), Hwang 
& Noh (2007), and Hwang et al (2012). 

We anticipate potentially wide applications of our ex- 
act and fully nonlinear perturbation formulation, not 
only in higher order perturbation theory, but also in 
the averaging, fitting and back-reaction approaches in 
theoretical cosmology (Ellis 1984; Ellis & Stoeger 1987; 
Clarkson et al 2011). 
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KOSEF funded by the Korean Government (MEST). 
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APPENDIX A. ADM (3 + 1) EQUATIONS REVIEW 

The ADM (Arnowitt-Deser-Misner) formulation (Arnowitt et al. 1962) is based on splitting the spacetime into the 
spatial and the temporal parts using a normal four- vector field n a . The metric is written as 

g m = -N 2 + N l N l: g Ql =N t , g,, /,,,.. g 00 = -N~ 2 , g° l = N~ 2 N\ <f J = h lJ - N~ 2 N l N J , (92) 

where N% is based on hij as the metric, and h 1 ^ is an inverse metric of hij; for meanings of the ADM variables, see 
Smarr & York (1978). The normal four- vector n a is introduced as 

n = —N, ni = 0, n° = TV" 1 , n l = -N^N 1 . (93) 
The fluid quantities arc defined as 



E = n a n b T a \ Ji = -n b T?, S l3 =T. l3 , S = h ij S i:j , Sy = Sy - (94) 



where Ji and Sij are based on hij. The extrinsic curvature is introduced as 



^ = ^(^+^-y. K = h*iK iS , K, s K,, ^h,,K. i or, i 



where is based on hij. A colon ':' denotes a covariant derivative based on hij. ^ h ^ k is the connection based on 
hij as the metric, r ^ = \h lt (hji,k + hki,j ~~ hjkj). The intrinsic curvatures are based on hij as the metric 

T}{h)i _ -p(h)i _-p(h)i -p(h)m-p(h)i _ -p(/i)m -p(h)i 

n jhl — 1 jt,k 1 jk,t <~ 1 jl L km 1 jk L Cm' 

Rf^R^\ kj , flW S A«J#\ R^^Rf-\ hij R^. (96) 
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A complete set of the ADM equations is the following (Bardeen 1980; Noh & Hwang 2004) 



RW = K l3 K %0 - -K 2 + l6irGE + 2A, (97) 



K\. } - ^K,t = &irGJ t , (98) 

A oA" 1 - K^N'N' 1 + A : %A _1 - Tt^Kij - -K 2 - AttG (E + S) + A = 0, (99) 

3 

K'jflN- 1 -K i j:k N k N- 1 + K^N'.^N- 1 -T? k N k . j N- 1 = KK 1 ^ - (n'^ - ~5}N :k k ^ A" 1 + R (h) ) - 8ttG^,(100) 

EflN- 1 - E^N'N- 1 - K (e+'S^J -^Kij+N' 2 (A 2 J*). . = 0, (101) 

J^oA" 1 - J^WN- 1 - JjW.iN- 1 - KJi + EN^N- 1 + Sj :j + SfNjN -1 = 0. (102) 

These are the ADM energy constraint, the ADM momentum constraint, trace of ADM propagation, tracefree ADM 
propagation, the ADM energy conservation, and the ADM momentum conservation equations, respectively. Equations 
(|97[) - (|102l) together with definition of K in equation (I95[) provide the fundamental perturbation equations presented 
in Equations (|2"Tj) - (l2~T|) . In the Appendix B we present details of steps useful for the derivation. 

APPENDIX B. DERIVATION OF FULLY NONLINEAR PERTURBATIONS 
Our metric convention is 

5oo = -a 2 (1 + 2a) , g i = -axi, 9ij = a 2 (1 + 2ip) % , (103) 

where Xi IS based on 5ij as the metric. We assume a is a function of conformal time (x° = rf) only, whereas a, Xi an d 
tp are general functions of space and time, but we do not assume these to be small in amplitudes. Justification of our 
metric convention is made in Section [5] The inverse metric is 

~oo _ __L l + 2<p 0l _ 1 xVa 

9 ..9/1 i r,._\/i , n„\ , ..fc.. 9 



a 2 (1 + 2^(1 + 2a) + X k Xk/ 'a 2 ' » a 2 (1 + 2^(1 + 2a) + X k Xk/a 2 ' 

qV = ( X l X J /a 2 \ , 1Q4 v 

9 a 2 (l + 2<p) \ (1 + 2^(1 + 2a) + X k Xk/a 2 J ' 1 ' 



From Equation (|92|) the ADM metric can be identified as 



thus 



and 



hij -a 2 (1 + 2^)%, /^ = -_L_^, (105) 



^r00 _ \ ~0i X ~ij _ 1 ( cij X X" 1 i I IOC I 

y aW 2 ' 9 oW 2 (l + 2^)' y a 2 (l + 2^) v a 2 AA 2 (1 + 2<p) • ' 



n, = 0, n = -aN, n l = „ *' - . , n° = -\- r (107) 



aW(l + 2^)' aA" 



The three-space connection and curvatures are 



r (h)i _ 
1 jk ~ 


1 


l + 2<p 


flW = - 




l + 2cp 


# j - 


1 


a 2 (l + ^ 



o f A i y y.fe \ s R (h) - 2 f o a M i ? ^ y '- fc 

1 + 2^ (l + 2^) 2 + (1 + 2^)V iJ ' a 2 (l + 2<^) 2 V ^ 1 + 2^ 



(108) 
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The extrinsic curvature gives 

„2 



(if + £ + 2H<p) % + ^ (xij + Xj,i) 



A 
K 



1 



1 



a 2 (l + 2y>) 



X fe ¥>,fc 



(Xi<P,j +Xj<P,i -X k V,k5i 3 ) 



Af(l + 2tp) 
1 

' ~a?Af(l + 2<p) 
1 



3 (ff + ^ + 2^) + ^X fe . fe nr 



-3F + k, 



5^- 



1 



7f = 

J * aW 2 (l + 2(^) 2 



1 ; , , 1 . 4 



1 + 2^ 

i , J 



2 

3 Sv 



2% lJ (XiJ + Xi,i) - gX 1 , ( X J j - y 



2(^ [2 



t:xV j (xij + xu) - \x l ,,xVj 



(l + 2(^) 2 

The fluid four- vector is 



X^W 3 <P,j + ^'xV.^j 



(109) 



AMI 



l + 2yj 1 + 2^ 



v 



1 



X 1 L v k Vk ~o _ 

a(l + 2tp) a 2 Af(l + 2ip)\ l + 2p' aA/" 1 



T 



v k v k 



(110) 



1 + 2^' 

where Uj is based on <5y as the metric. For Uj = we have u a = n a . The energy- momentum tensor of an ideal fluid is 

1 



H+P 



1 + 2^ 

= a 2 [(1 + 2<p) pSij + (fi + p) ViVj] 
The ADM fluid quantities become 



v k Vk 
l + 2ip 



rpK} 
i 1 i 



Af 



v k Vk 
l + 2<f 1 



(111) 



V k Vk 



E = (i + ((j,+p) , J i =a((Jt + p)v i \l + 



l + 2<p 



V k Vk 

l + 2<p 



a(l + 2^) 



v\l + 



V k Vk 

1 + 2^' 



1 + 2lp 



, S = 3p- 



M + P 
l + 2<p 



v Vk, Si = 



l + 2<p 



(112) 



Notice that for Vi — (ignoring the vector-type perturbation = 0, and taking the comoving gauge v = 0) the 
energy-momentum tensor and the ADM fluid quantities are simplified a lot with 



T ° = -ju, T° = 0, T l3 = a 2 (1 + 2tp) pSij , # = Jl, J; = 0, S = 3p, 5. ■ = 0. 



(113) 



Using the above quantities, from A in Equations ([55)1 and (IWl) - (|102p we can derive Equations ([2"TT) - (|27p . respectively. 



APPENDIX C. COVARIANT FORMULATION 



Here we present the covariant kinematic quantities based on the normalized fluid four- vector u a . We have (Ehlcrs 
1993; Ellis 1971, 1973) 



h a h b u c -d = h[ a h b ]U C ;d + h( a h b yu c; d = uj ab + 9 ab = u a . b + a a u b , 
&ab = Qab - ~0h ab , 9 = u a . a , a a = u a = u a;b u b , 



(114) 



where h ab = g a b + u a u b is the projection tensor with h ab u b = and h a a = 3; t( ab ) = \{t a b + ha) and t[ ab ] = \ {t ab ~t ba ); 
we have a c u c = and a ab u b = = Lu ab u b . An overdot with tilde ' indicates a covariant derivative along u a . The 
quantities 9, a c , aij, and u>ij are the expansion scalar, the acceleration vector, the shear tensor, and the rotation 
(vorticity) tensor, respectively, of the u a flow. 
Using Equations ([92]) . (|105[) . (|108[) and (|110[) we can show [the spacetime connection T bc in terms of the ADM 
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notations is presented in Equation (6) of Noh & Hwang (2004)] 



= (3H-k) Wl 



v k v k 



(A/V),. 



+ 



1 



9o 



X' 



1 + 2<^ a7V(l + 2<p) o(l + 2^) 2 aA/" V a (! + 2( f) 



Vih/1 



v k Vk 
1 + 2^' 



\ + 2<pJ N l + 2ip (l + 2ip) 2 aAf\ 1 + 2^ 
= au[i,j] - a«[i{ fl + 
1 



/ \ , ^,fcX fc , / x fc 



2^5 



l + 2ipj AT l + 2ip (l + 2^) 2 



1 



aM V 1 + 2<£ 



, «j],fcX fc , / x fc 



l + 2(p 



— / V L Vl 

~^V 1 + TT2^ + a ^ ) 

^j),fc^ fc vkv k<P,j) 1 

"l + 2v? ~ (l + 2(p) 2 + ^A7V X T l + 2ip 



v*-vg 



/ \ , ^),fcX fc / x fc 



i 



3 
1 a 



a(l + 2<p) 

[(1 + 2</j) <5jj + ViVj 



, , (3fr-K)^/i 



v k v k 



l + 2ip 



l + 2(p 



X 



a(l + 2ip) 



V fc 4 1 



l + 2<£ 



(115) 
(116) 



(117) 



(118) 



where do = J^. Taking a normal frame with w a = n a , we have i?j = 0, thus 



'(n) _ 



= -K = 3H 



-in) -A^i 



-0, <7y 



= -if,: 



(119) 



The same is valid for an irrotational fluid in the comoving gauge so that we have Uj = 0. Notice that even for vanishing 

vector-type perturbation with = = ^ , we have ^ to nonlinear order in many other gauge conditions; in 
the comoving gauge we have Uj = 0, thus u>ij = 0, see Section [7] 
Part of the covariant equations are (Ehlers 1961; Hawking 1966; Ellis 1971, 1973) 



-0 2 - a a . a + a ab a ab - u ab u ab + 4irG (Ji + 3p) - A = 0, 



~ __h b a p,b 



(120) 
(121) 
(122) 



These are the Raychaudhury equation, and the covariant energy and the covariant momentum conservation equations, 
respectively; the latter two equations are presented in the energy- frame with vanishing energy flux q a = (Ellis 1971, 
1973; Hwang & Vishniac 1990); the fluid four vector u a in the energy frame (q a = 0) coincides with the normal frame 
four-vector n a for = (the irrotational and the comoving gauge). By combining Equations (|120[) - (|122[) we have 
(Jackson 1972, 1993; Hwang & Vishniac 1990) 



J^— = 4irG (/7 + 3p) - A + U 2 + a ab a ab - u ab u ab - a c a c + h ab ( , 
fj, + p 3 \H+Pj. b 



(123) 



Using the covariant kinematic quantities based on u a presented above we can derive nonlinear perturbation equations 
in gauge-ready forms. Equation (|120[) leads to a complicated combination; in the normal frame it gives Equation (|24|) 
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A7 1 



l 



v k v k 



d 



l + 2(p \dt a 2 (l + 2ip) 

(A/V),, v i (f 



V, )Jl + (ji+p)(3H-K)\ 1 + 



1 



8 



X' 



A7 1 



aAf(l + 2ip) a(l + 2ip) 2 Af \dt a 2 {l + 2ip) 
V 



v k Vk 
l + 2ip 

V, 



a(l + 2<p) 



1 



9 



1 + 2tp \dt a 2 (l + 2tp) 



1 



1 



V k Vk \ AT, 



\ + 2<p) J\f 



aAf^ 



1 



l + 2(^ V 1 + 2 V 3 



P. 



Xj 



N 1 



1 



V k Vk 

l + 2ip 



d 



(124) 



X'' 



1 + 2^ V^t ' a 2 (l + 2ip) 
v k v k (p,i 1 



ViV k p,k 



Vi,kV 

l + 2(p^ {I + 2(p) 2 l + 2(p JI + p 



(125) 



These are alternative forms of Equations (|26|) and ([27]) , respectively. The above two equations follow from the covariant 
energy and momentum conservation equations based on projecting along u a and h ab as 



= T b u a 



_ b . a = + = T^ b h c a = (ji + p)Z a + h b a p, b . (126) 

Whereas the ADM energy and momentum conservation in Equations (|2"B)) and (|2"T|) are based on projecting along n a 

~{n) 

and h , = g ab + n a n b . We have the following relations 



1 



N 



1 



1 u 



= Tl^ = -T b :b —T b b = -- (no - N%) T* ]b u° - - ( - + AT ) T^ b h 



N 



[J — 1 a;b n i — 1 i;b — 1 a;b n i u i J a;b u • 



N 

^b 



N \ u 



-ib ua 



These lead to Equations (fToTj) and (fT02)) . 
A covariant spatial gradient variable 



A a = l:h b a n,b, 



(127) 



(128) 



was introduced as a covariant and gauge-invariant variable in the linear perturbation theory (Ellis & Bruni 1989; 
Woszczyna & Kulak 1989). To the fully nonlinear order, we have 



A, = i 



Vi Jrt 



v k v k 



d 



l + 2ip \dt a 2 (l + 2^) 



2tp 



(129) 



Using the gauge transformation properties to the second order in Hwang et al (2012) we can show that A, is gauge- 
invariant only to the linear order. To the linear order, we have Aj = S V i + wj which is a sum of the spatial gradient of 
S in the comoving gauge and the vector-type perturbation. In the comoving gauge without vector-type perturbation 

(thus Vi — 0) we have Aj„ = S v ,i/(1 + S v ) which is related to the density perturbation in the comoving gauge. 

APPENDIX D. MULTIPLE COMPONENTS OF IDEAL FLUID 

Here we consider the case of multiple components of ideal fluid. Even in the presence of many fluids (with vanishing 
fluxes and anisotropic stresses) all our equations are valid with the fluid quantities considered as the collective ones. 
In the presence of N fluids we have 



Tab - T(Jj ab , 



(130) 



with the fluid quantities of collective and individual components introduced as 

T ab = Jm a u b + p(g ab + u a u b ) , T {I)ab = Ji(i)U (I)a u {I ) b + p (I) (g ab + u (I)a u {I)b ) , 
where indices /, J, • • • = 1, 2 ... N identify the fluid component. From Equation (|130l) we have 



M - 3 P = X! - 3p ( j)) , u a 



u (j) u t 



(131) 

(132) 
(133) 
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Now we introduce the normalized {u c ^u^iy = — 1) fluid four- vector of individual component as 

= av(i)i- (134) 

the rest of W(/) a are the same as in Equation (|110|) with replaced by ivm. Similarly, for the energy- momentum 
tensor, the ADM fluid quantities, and the covariant kinematic quantities of the individual fluid component, we can 

replace Jl, p and Vi to Jim, pm and vmi, respectively, in Equations (jlll|) - ()118j) with T a b, E, 9, etc., replaced by Tm a b, 

E(i), etc., respectively. For the ADM fluid quantities we have 



(J)j> 



(135) 



but these simple relations do not hold for the covariant kinematic quantities 9, a.;, u>ij and cr^-. From Equation (|110p 
for the collective component and for the /-component, we have 



u (I) u c — 



v k {I) v k 



1 + 



2ip \ 



l + 2<p 



v'-vi 
l + 2<p 



Equation (|133|) gives 



(136) 



(137) 



Equations (j 132[) and (|137|) with Equation (|136[) provide relations between the collective and indivdual fluid quantities. 
Using these relations our fully nonlinear and exact perturbation equations remain valid even in the multiple component 
fluid case. 

Now, we have to provide the equations followed by the individual fluid. The energy and the momentum conservation 
equations for individual component are the ones we need. The energy and momentum conservation equations follow 
from f Q b . b = 0, thus 



T (I)a;b = I (I)a, /] I(J)a — 0; 



(138) 



where I(i) a is the interaction terms among fluids. In the ADM and the covariant formulations, the energy and the 
momentum conservation equations are presented in Equations (47) and (48) of Noh & Hwang (2004) , and Equations 
(31) and (32) of Hwang & Noh (2007), respectively. In ideal fluids, these are 



E^N- 1 - Em^N- 1 - K (e {i) + ifyj) - S^K,^ + N~ 2 (n 2 ^), = ~ (lm - h^ N '' 
Jd^oN- 1 - Jmi-.jWN- 1 - JmjW.jN- 1 - KJ {I)l + EmNiN' 1 + S (I j i:j + S^NjN- 1 = 

hn + (hi) +P(i)) = -u a (i)T(i)a, 

h (I)aP(I)-b _ h (I)a I (I)b 



(139) 
(140) 
(141) 

(142) 



Compare these with Equations (|101j) . (I102p . (|12ip , and (I122[) , respectively, for collective ones. From Equations (jl39p ~ 
(|142p we can derive the energy and the momentum conservation equations and the alternative forms. The results 
are the same as Equations (|26|). (|27)l . f!124[) and (|125[) with all fluid quantities Ji, p and V{ replaced by Jim, P(i) and 
v/ju, respectively, and add the following contributions from interactions among fluids to the right-hand-sides of the 
equations 



1 (~ X 1 



a(l + 2(p)' 



■a(l + 2 ¥J ) U « 7 W-^V 1 
1 



W) fc J x * r 



m +p(i) 



a(l + 2ip)' 
v (i) v (i)k 



T+Tpv 1 ™ + J! V 1 + TW l /(/)0 + 



l (i)j 



(143) 
(144) 

(145) 
(146) 
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respectively; we have introduced I<i)i = where the spatial index of is based on 7^; Inn is the perturbed 

order quantity whereas I(i)o can include background order quantity as /(/)o = -^(7)0 + SI(i)o- These equations together 
with Equations (|132j) and ()137j) complete the additional equations we need in the multiple fluid system. 
The vector variables V(i)i and Itni can t>e decomposed to the scalar- and vector-type perturbations as 

V(i)i = + vffi v = 5/(j),j + 6Iff v (147) 

As in the single component case, to the linear order the scalar-type perturbation 6n), Spii), &I{i)q and 51^) 
depend on the temporal gauge transformation whereas the vector-type perturbations are gauge invariant (Hwang 
1991). In addition to the fundamental gauge conditions in Equation (|18[) . in the multi-component case, we have the 
following gauge conditions available 

I— component — comoving gauge : Vn) = Qj 

uniform — I — component — density gauge : 6^=0, (148) 

to the fully nonlinear order. 



REFERENCES 



Arnowitt, R., Deser, S., & Misner, C. W. 1962, Gravitation: an 
introduction to current research, edited by Witten, L. (Wiley, 
New York) 

Bardeen, J. M. 1980, Phys. Rev. D, 22, 1882 

Bardeen, J. M. 1988, Particle Physics and Cosmology, edited by 

Fang, L., & Zee, A. (Gordon and Breach, London) 
Clarkson, C, Ellis, G., Larena, J., & Umeh, O. 2011, Rep. Prog. 

Phys., 74, 112901 
Ehlers, J. 1961, Proceedings of the mathematical-natural science 

of the Mainz academy of science and literature, Nr. 11, 792; 

English translation in 1993, Gen. Rel. Grav., 25, 1225 
Ellis, G. F. R. 1971, General relativity and cosmology, 

Proceedings of the international summer school of physics 

Enrico Fermi course 47, edited by Sachs, R. K. (Academic 

Press, New York) 
Ellis, G. F. R. 1973, Cargese Lectures in Physics, edited by 

Schatzmann, E. (Gordon and Breach, New York) 
Ellis, G. F. R. 1984, General Relativity and Gravitation, edited 

by Bertotti, B. et al. (Reidel, Dordrecht) 
Ellis, G. F. R., & Bruni, M. 1989, Phys. Rev. D, 40, 1804 
Ellis, G. F. R., & Stoeger, W. 1987, Class. Quant. Grav., 4, 1697 
Field, G. B., k, Shepley, L. C. 1968, Astrophys. Space Sci., 1, 309 
Friedmann, A. A. 1922, Zeitschrift fur Physik, 10, 377; English 

translation in 1999, Gen. Rel. Grav., 31, 1991 
Harrison, E. R. 1967, Rev. Mod. Phys., 39, 862 
Hawking, S. W. 1966, ApJ, 145, 544 

Hawking, S. W., & Ellis, G. F. R. 1973, The large scale structure 

of space-time (Cambridge University Press, Cambridge) 
Hwang, J. 1991, ApJ, 375, 443 
Hwang, J. 1994, ApJ, 427, 542 

Hwang, J., & Noh, H. 1999, Gen. Rel. Grav., 31, 1131 
Hwang, J., & Noh, H. 2001, Phys. Rev. D, 65, 023512 



Hwang, J., & Noh, H. 2005, Phys. Rev. D, 71, 063536 

Hwang, J., & Noh, H. 2006, Phys. Rev. D, 73, 044021 

Hwang, J., & Noh, H. 2007, Phys. Rev. D, 76, 103527 

Hwang, J., Noh, H., & Gong, J. 2012, ApJ, 752, 50 

Hwang, J., Noh, H., & Park, C-G. 2010, Phys. Rev. D, 82, 043525 

Hwang, J., & Vishniac, E. T. 1990, ApJ, 353, 1 

Jackson, J. C. 1972, Proc. Roy. Soc. London A, 328, 561 

Jackson, J. C. 1993, MNRAS, 264, 729 

Jeong, D., Gong, J., Noh H., & Hwang, J. 2011, ApJ, 727, 1 
Kodama, H., & Sasaki, M. 1984, Prog. Theor. Phys. Suppl., 78, 1 
Landau, L. D., & Lifshitz, E. M. 1975, The Classical Theory of 

Fields (Oxford, Pergamon) 
Lifshitz, E. M. 1946, J. Phys. (USSR), 10, 116 
Lukash, V. N. 1980a, Sov. Phys. JETP Lett., 31, 596 
Lukash, V. N. 1980b, Sov. Phys. JETP, 52, 807, 
Lyth, D. H., Malik, K. A., & Sasaki, M. 2005, JCAP, 05, 004 
Ma, C, & Bertschinger, E. 1995, ApJ, 455, 7 
Mukhanov, V. F. 1988, Sov. Phys. JETP, 68, 1297 
Mukhanov, V. F., Feldman, H. A., & Brandenberger, R. H. 1992, 

Phys. Rep., 215, 203 
Nariai, H. 1969, Prog. Theor. Phys., 41, 686 
Noh, H., &; Hwang, J. 2004, Phys. Rev. D, 69, 104011 
Peebles, P. J. E. 1980, The Large-Scale Structure of the Universe 

(Princeton University Press, Princeton) 
Sachs, R. K., & Wolfe, A. M. 1967, ApJ, 147, 73 
Sasaki, M. 1986, Prog. Theor. Phys., 76, 1036 
Smarr, L., & York, J. W. 1978, Phys. Rev. D, 17, 2529 
Vishniac, E. T. 1983, MNRAS, 203, 345 

Woszczyna, A., & Kulak, A. 1989, Class. Quantum Grav., 6, 1665 
York, J. W. 1973, J. Math. Phys., 14, 456 
Zel'dovich Ya. B., & Novikov, I. D. 1983, Relativistic 

astrophysics, Vol 2, The structure and evolution of the universe 

(Univ. Chicago Press, Chicago) 



